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Abstract 

We consider the time dependent probability distribution of a coarse grained observ- 
able Y whose evolution is governed by a discrete time map. If the map is mixing, 
the time dependent one-step transition probabilities converge in the long time limit 
to yield an ergodic stochastic matrix. The stationary distribution of this matrix is 
identical to the asymptotic distribution of Y under the exact dynamics. The nth 
time iterate of the baker map is explicitly computed and used to compare the time 
evolution of the occupation probabilities with those of the approximating Markov 
chain. The convergence is found to be at least exponentially fast for all rectangular 
partitions with Lebesgue measure. In particular, uniform rectangles form a Markov 
partition for which we find exact agreement. 
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1 Introduction 



The technique of coarse graining has often been used to obtain irreversible be- 
havior from systems which, though they may be reversible, are "sufficiently" 
chaotic [7,12,15]. The procedure is to describe the time evolution of the sys- 
tem in terms of macroscopic variables whose preimages form the cells of the 
partition. More abstractly, coarse graining may be described as a contraction 
of the exact dynamics onto a corresponding symbolic dynamics [1] . One of the 
goals of this approach is be able to describe the macroscopic dynamics in a 
self-contained form, in other words, to obtain macroscopic laws which do not 
refer to the underlying microstates. This may be possible if the macroscopic 
dynamics satisfy a master equation [13]. 
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It has long been realized, however, that the evolution of observables under 
reversible microdynamics are generally non-Markovian due to memory effects 
[12]. While it may be possible to construct a partition for which the resulting 
symbolic dynamics are exactly Markovian [15,11], such partitions are excep- 
tional and in general do not correspond to the macroscopic variables one is 
interested in. Assumptions of molecular chaos have in the past been used to 
argue that the Markov assumption is well approximated for complex macro- 
scopic systems [12]; however, it is unclear whether such assumptions would be 
consistent with the underlying microscopic dynamics. What are needed are 
sufficient conditions on the dynamics, as well as the macroscopic observables, 
to ensure that the Markov assumption is well approximated. This has been 
done, e.g., for quantum systems weakly interacting with a heat bath [5], but 
it is not clear how those results may apply in more general settings. Since 
the master equation is characterized by only short-term memory effects, it is 
natural to suppose that decaying correlations should at least be a necessary 
condition. Thus, one is led to consider systems which are mixing. 

We begin in Section 2 with a discussion of the general problem for discrete time 
maps. An equation for the time dependent occupation probabilities, i.e. the 
generalized master equation, is derived. (Since it involves one-step transition 
probabilities that are in general dependent upon both the initial cell occupied 
and the iteration time, the generalized master equation is not generally a true 
master equation.) In Section 3 we consider the asymptotic behavior of these 
transition probabilities for long times. For mixing maps it is shown that the 
one-step transition probabilities converge to an ergodic stochastic matrix. The 
corresponding Markov chain is found to converge to the same stationary distri- 
bution as that of the exact process. The baker map is considered in Section 4 
as an exactly solvable system which nevertheless contains many of the impor- 
tant features of more physical systems. An explicit expression is obtained for 
the nth iterate of a point under this map, thus allowing us to compute exact 
time dependent occupation probabilities. For rectangular cells with Lebesgue 
measure it is shown that these probabilities converge at least exponentially 
fast. In Section 5 we compare the occupation and transition probabilities de- 
termined by the baker map with their corresponding approximations using the 
approximating Markov chain. In particular it is shown that for the baker map 
any uniform rectangular partition is a Markov partition for which the Markov 
approximation is exact. The main conclusions are summarized and discussed 
in Section 6. 



2 Generalized Master Equation 



Let {X,A,P) be a probability space and let be a measurable map on X. 
We assume that P has a density p with respect to some finite measure //, i.e. 



P ^ fi, and that is /i-measure preserving, i.e. /i o 0~^ = /i. For xq G X and 
n G Z, 0"'(a;o) is the n^^ iterate of Xq under 0. (If is invertible, then 0" = 
(0~^)l"' for n < 0; otherwise, 0" is undefined.) The preimage of a measurable 
set A under 0" will be denoted by (f)~"'{A) and represents the set of points 
that map into A after n time steps. Thus, P[0""'(A)] is the probability that 
an iterate will be in A at time n. 

We wish to consider the evolution of a phase function Y of the form 

Yix)^J2y^'^cAx): (1) 



where / is a finite index set, {i/jjje/ are distinct real numbers, and Ic^ is the 
indicator (characteristic) function on Ci G A. Without loss of generality, we 
may assume {Ci}i^i is a disjoint partition of X. We further suppose that the 
partition is nontrivial in the sense that /u[Cj] > for all ? G /. Let y„ = F 00" 
and define the transition probabilities T/"^ and T^™'" as follows: 



<n) 



T,7^P[0-"(C,)|a], (2) 

T^f = p[r'\c,) 1 0-(c,) n a] , (3) 

for P[Ci] > and P[0~'"(Cfc)nCj] > 0, respectively, and zero otherwise. We see 
that Tjj is the probability that F„ = yj given that Yq = Hi, and TJ^J^j is the 

probability that F„ = yj given that Ym = Vk and Fq = Hi- In general, T^^j 
will depend upon i, so that {Yn}nei is generally a non-Markovian random 
process. 

Using Eqs. (2) and (3) we obtain the generalized Chapman- Kolmogorov equa- 
tion for this process: 

rp{n+m) ■^^ rp{n)rp{n,n+m) r ,\ 

k 



Let W^l^j = T^^'"^ denote the one-step transition probability at time n con- 
ditioned on Ci. We see that W^J^j is the probability that an iterate in cell Ck 
at time n will make a transition to cell Cj at the next time step, given that it 
started in cell Cj. Using Eq. (4) it can be shown by induction over n that 

n 

T!f = T.---T. liw^Z-l^^ (5) 

fcl fcn-l in=l 



where ko = i and kn = j- Let Pj{n) = P[0 "(C^)] denote the occupation 



probability of cell Cj at time n and note that 

p,{n)=J:P^WTif . (6) 

i 

This is the generalized master equation for the exact occupation probabilities. 
In general, W^J^- depends on both i and n, so Eq. (6) contains memory effects, 
i.e. the process is not Markovian. Nevertheless we might expect quasi Marko- 
vian behavior for maps with decaying correlations. This suggests that mixing 
may be a key property to consider. 

3 Mixing and Approximating Markov Chains 



The map defined in Sect. 2 will be mixing with respect to n [10] if for all 
measurable sets A and B, 

lim /i[0-"(A) nB] = P,[A]ii[B] , (7) 



where P* = -jx^fJ'- From Eq. (7) it follows [10, p. 72] that if P -< /i and 
P[P] > then 

limP[0-"(A)|P] = P4A]. (8) 



Using Eq. (8) we find that for P[Ci] > 0, 

isi^5^ = iL^^r"(^.) I c^] = p.[c,] , (9) 



the limit being zero otherwise. Note that from Eqs. (6) and (9) it follows that 
Pj{n) -^ P*[Cj] as n ^ oo. Thus, mixing implies an asymptotic convergence 
of solutions to Eq. (6) for any initial P ~< n. From Eq. (8) we also note that 
if P[Ci] > then P[0-"(Cj) n C^] > for all n sufficiently large. (Recall that 
P*[Cj] > for a nontrivial partition.) Thus for P[Ci] > we similarly have 

hm T/-r-) = hm ^^'^"^""^^.'^•\'^^',^'^^^ = P.[r-{C,) I C] . 

n^oo *'^J n^oo P[0-"(Cfc) | Cj] ^ ^ 



The one-step transition probabilities therefore have the following asymptotic 
form, independent of k: 

lim l^g. = P*[<f^'\C,) I a] ^ W,, . (10) 



Since W is clearly a stochastic matrix, it must correspond to some Markov 
chain {Yn}nez on the state space {yi}i£i. We shall refer to {Yn}nez as the 
approximating Markov chain for the exact process {Yn}n(^i- Let qi(n) denote 
the probability that Yn = Hi for i G /. In contrast to Eq. (6), we have simply 
[13, p. 95] 

g,W = E^^(o)(w^")M- (11) 



The condition T/^ = {W"')i^j is thus seen to be the necessary and sufficient 
condition for exact agreement between {Yn} and {Yn}, assuming j9j(0) = (^^(O) 
[11]. The following theorem establishes that, even without this equality, the 
two processes converge to the same stationary distribution. 

Theorem 1 If (p is mixing with respect to a finite invariant measure /x and 
{Ci}i(zj is a nontrivial finite partition of X , then any approximating Markov 
chain is ergodic. 



PROOF. A Markov chain is ergodic if it is irreducible, aperiodic, and pos- 
sesses a stationary distribution [16, p. 109]. Let i,j G / be fixed and take 
P = P* in Eq. (5). Given that is mixing, we know Tj" -^ P^[Cj\. Since 

IJi.\Cj\ > by assumption, P*[Cj] > and so T/" > for all n > Nij, where 
Ni j is some positive integer. Thus, 



ki,...,k„-ieK„ m=l 
1 " 

< - E n P*[<f^'"'iCkJ n (/^-^""-'HCk^^J n a] , 

P" ki,...,k„^lGK„m=l 

where ki = i, kn = j, and 



Kn = \ku..., kn-i G / : n n[0-('"-'n^fc^-i) n a] > , (12) 

I m=l J 

Pn = ^ min |nP40-(™-i)(C,_Jnc4 • (13) 

We are assured that Kn is nonempty since, were it empty, Tj " would be zero. 

Since we have taken n to be sufficiently large so that T^J > 0, it follows that 
Kn is indeed nonempty, and so p„ > 0. Multiplying both sides of the inequahty 



by pn we conclude that 



kl,...,kn-l£Knm=l 



Now, since Wij > P*[0 ^{Cj) fl Cj] and P* is invariant, 



ki fc„_i m=l 



ki fcn-l "1=1 

=E--- E n ^*r'"(c.jn0-(-^)(c._j] 

n 

> E n P*[r"\C,J n 0-("-^)(C^fc™-J n a] . 

ki,...,kn-ieKn m=l 

It then follows that {W^)ij > for n > Nij. Since i and j were arbitrary, this 
proves irreducibility. 

A Markov chain is aperiodic if it has period d = 1. If we suppose that d > 1, 
however, then for 1 < r < rf we would have that {W"^'^~^^)i^i = for all m G N. 
We have just shown, however, that (W^")i,i > for all n sufficiently large. 
Since d cannot be smaller than 1, we conclude that d = 1. 

Finally, it may be readily verified that qi{0) = P^[Ci] for i G / is a stationary 
state of W. □ 



An ergodic Markov chain will converge to a unique stationary distribution. 
We thus conclude that for all i,j 

\ini{Wn., = P.[C,]. (14) 



The limit is both geometric in rate and uniform in i,j [4, p. 124]. We conclude 
that both qj{n) and Pj{n) approach the same asymptotic value of P^,[Cj]. The 
asymptotic distribution will be uniform whenever the partition is uniform, 
i.e. when fi[Ci] = fJ^lCj] for all i,j G /, and the following theorem gives the 
necessary and sufficient conditions for this to be true. 

Theorem 2 The partition is uniform if and only if W is doubly stochastic. 



PROOF. If the partition is uniform then 



Now suppose W is doubly stochastic and note that for any n G N, W^" will 
also be doubly stochastic. Since the partition is finite, 

1 = hm E(w^")m = E 1™ (w^")„- = E^4Q] • 



If there are M partition cells, then ^[Cj\ = j^fi[X] for all j E I. □ 



4 The Baker Map 



We now consider the above results in light of a particular example, the baker 
map, which is defined on X = [0, 1)^ as follows [10]: 

(l){x,y) = [2x, -j l[o,i/2)(a:) + (^2x - 1, - + -j l[i/2,i)(a;) . (15) 



The baker map is a Lebesgue measure preserving map of the Bernoulli type 
and hence is mixing [2]. It is also time reversal invariant in the sense that 
(j) o R = R o (j)~^ for the transformation R{x,y) = {y,x). Thus, it possesses 
all the salient features of a closed Hamiltonian system with time reversal 
invariance while permitting a simple and exact analysis. An example of an 
real-world optical system which produces a baker transformation is described 
in [8]. 

We wish to explicity compute the probability Pj{n) defined in Section 2. To 
do so we will need to calculate probabilities using the given initial probability 
density. Now, given a probability density p, the density after n time steps 
is given by U^p, where U is the Perron-Frobenius operator. Although the 
generalized spectral decomposition of U (and hence U"') for the baker map has 
been constructed [9], it cannot generally be used for computing probabilities. 
However, since is invertible and area preserving, f/"p = p o 0~" [10, p. 47]. 
Furthermore, since 0~" = i? o 0" o i?, to compute f/" for n G Z it suffices to 
determine 0" for n G N. 

Theorem 3 Let denote the baker map and let U be the corresponding 



Perron- Frobenius operator. For any integer n ^N and point {x,y) G [0, 1)^, 



2" 



r{x, 



,y)= J2 (2"x-m + l, |^ + a(m)jlr_i „.)(a;) , (16) 

m=l ^ ^ / 



where 



,0, if m = 1, 

a{m) = { (17) 

a(m - 2"-i) + 1/2", zf m e {2"-^ + 1, . . . , 2"} . 



PROOF. The proof is by induction on n. For n = 1 Eq. (16) clearly holds. 

By Eq. (15) we have, for 0"(x, y) = 0(0"~^(x, y)). 



on— 1 

0"(x 



,1/)=^ (2"a;-m + l, — + a(m)j lj^__^^)(2x) l[o,i/2)(a;) 

2""i / 1 \ 

+ V (2"x - 2"-i -m + 1, ^ + — + a(m) ) 

m=l ^ Z, Z, / 

X Irm-i ^2_'\(2a;-l) lri/2i)(a;) 

In the second sum we make the substitution m = m' — 2"^^ and note that 
2^ + a{m) = an{rn'), by Eq. (17). Finally, we note the following identities: 

l[irL^,_irL^)(2a:) l[o,i/2)(a;) = l[i|^, ™ )(a:), 
since m < 2"'^^ is equivalent to m/2" < |, and 



lr m-1 m \(2x — 1) lri/2 i)(a;) = l^ m'-i m'Jx), 

[2^ — 1 ' 2"^^ / L / 1 / L 2" ' 2" ' 

since m' > 2"^^ + 1 is equivalent to (m' — l)/2"' > |. Thus, 



0"(x,i/)=5: f2"a;-m + l, |^ + a(m)j ![_,„.) (a:) 

+ V f2"a;-m' + l, |- + a(m')Vr'-'-i -' 

m'=2"-i+l ^ ^ / L 2 2 

which is identical to Eq. (16). □ 



Note that the sum in Eq. (16) has only one nonzero term, namely, m 
mn{x) = [2"a;J + 1. Thus it may be written more compactly as 

r{x,y)= (2"a;modl, ^ + a(m„(x))) . (1 



Using Eq. (16) and the fact that f/"p = p o (J? o 0" o J?), we find that for a 
rectangular region C = [a,b) x [c, d), 

^r"(C^)]= T. I I Piu,v)dudv, (19) 



where 7„,i(c) = max{2"c, / - 1} - / + 1 and 6n,i{d) = min{2"rf, /} - / + 1. (We 
define m„(l) = 2".) Note that 'yn,i{c) = 2"c mod 1 for / = 'm„(c) and is zero 
otherwise, while Sn,i{d) = 2"'d mod 1 for / = rrinid) and is one otherwise. 

Let us now consider an initial density of the form 

p{x,y) = EP^(0)//iL[a] lcAx,y), (20) 



where the pj(0)'s are given constants and fiL[Ci] is the Lebesgue measure of 
cell Cj. (Of course, such a density may be obtained by pre-coarse graining an 
arbitrary p, in which case Pi{0) = Ju_ p{x, y) dx dy.) If the cells are rectangular 
with Ci = [ai,bi) x [ci,di) then 



where 



fn{a) = ^{mm{bi, a + 6,/2"} - maxja^, a + a,/2"}) , (22) 

gn,i = ^{mm{di,6n,i{dj)} - max{Q,7„,;(cj)}) . (23) 

(The function {} is defined such that 'd{x) = x for x > and 'd{x) = 
otherwise.) The occupation probability Pj{n) is now found using Eq. (6). 

As an example, consider p = — hrrlcy For n < \og2{l/ dj) we have 

, , 1 ^(6,- - 2"a,)^(rf,- - 2"c,) 1 , ,, 

P.» = ^ \u ; M '. ' < ^ • (24) 



2" {bj -aj){dj -Cj) ~ 2' 



9 



If Oj = Cj = 0, we have an initially exponential decay of Pj{n) before it 
converges to the asymptotic value fihlCj] = bjdj. 

This example suggests a general convergence rate which is at least expo- 
nentially fast. Let us consider the rate of convergence of T^j. Assuming 
n > log2[(6j — aj)/{hi — ai)] and n ^ \og2{dj — Cj)~\ we have 



id-i 



^^L[Ci] T>y = fn[ailc)] 9n,l, + f[aild)] 9n,la + J2 /n["(0] 

l=lc+l 



where Ic = mn{cj) and l^ = mn{dj). The sequence a{l),a{2), ... is dense in 
[0, 1), taking values which are uniformly distributed throughout the interval. 
Thus, we may approximate the sum by an integral for large n. 



J2 fuHl)] ~ 2"(rf, - Ci){d, - cj) J Ua) da , 



l=lc + l 



where we note that 



/n(tt) 



a - Oj + 6j/2", if aj - 2^ < a < aj - ^, 
{hj - aj)lT^ a ai- ^ <a <bi- ^, 
hi — aj/2" — a ii hi — i^ < a < hi — ^, 
, otherwise . 



2. 
2"' 



(25) 



Upon integrating we find that 



p(n) 



flL[Cj 



fn[a{lc)] gn,h + fn[a{ld)\gn^ ^ 2{hj - ttj) 



iJ'dCi 



2^fiL[Ci_ 



Thus we have convergence which is at least exponentially fast. (Recall that 
convergence for the approximating Markov chain is always exponential.) 



5 Comparison of Exact and Approximating Processes 



We have found that for mixing maps both Pj{n) and qj{n) converge to the 
same limit, P^:[Cj], as n — *> oo. If furthermore Pj{0) = qj{0) for all j G /, then 
we expect only transient deviations oi pj{n) from its Markovian counterpart 



10 



qj{n). The difference between these two probabihties may be measured by an 
average over all cells of the ^°° norms: 



1 M 



^p-TjY. s^P \Pj ("-) - ^J (^) I ' (26) 



M — neN 



J= 



where M is the total number of cells. Note that < Dp < 1, the value of 
Dp depending upon both the initial measure P and the partition chosen. In 
particular, if P = P* then Pj{n) = P*[Cj] = qj{n) and Dp = 0. 

We now return to the baker map and consider a partition of two cells: Ci = 

[0, s) X [0, 1) and C2 = [s, 1) x [0, 1). For s < |, the transition matrix given by 
Eq. (10) is 

2{l-s) 2(l-s) 



The transition matrix for s > | is obtained by interchanging s with 1 — s and 
taking the transpose. 

In Fig. 1 we display the average deviation Dp for < s < 1. The asymmetry 
in the graph is due to the fact that the initial measure is on Ci; a mirror image 
of it is obtained if C2 is used as the starting cell. For s near 0, ^, or 1 we note 
that Dp goes to zero, implying exact agreement between Pj{n) and qj{n) for 
all j and n. This can be understood as a consequence of the following theorem: 



Theorem 4 Let P,,, = [^, ^] x [^, ^] /or ^ e {1, . . . , M}, j G {1, . . . , iV} 
and M,N eN. Then {Rij} is a Markov partition for the baker map. 

A proof of this theorem is given in the appendix. Since the initial measure 
was taken to be uniform with respect to the maximum entropy measure hl, 
it follows that {Yn} is a true Markov chain [3]. 



6 Discussion 



We have shown that for mixing maps the time dependent occupation and one- 
step transition probabilities, Pj{n) and WJ^^j respectively, converge to values 
which are independent of both the index k and the initial probability den- 
sity p. The resulting asymptotic transition probabilities Wij form an ergodic 



11 



Dr, 




Fig. 1. The above plot shows the values of Dp for a partition of two vertical cells of 
width s and 1 — s. Note that Dp goes to zero near s = 0, 2, 1- 

stochastic matrix whose unique asymptotic distribution is identical to that of 
the exact, non-Markovian, process. 

To examine the detailed evolution oi pj{n) we considered the baker map 0. An 
explicit expression for 0"(a;, y) was found with which we were able to determine 
explicit formulas for Pj{n) in the case of an initial density which is uniform 
on rectangular cells. It was shown that in this case the occupation probability 
Pj{n) converges at least exponentially fast, i.e. no slower than qj{n). 

Finally, we studied the average deviation between Pj{n) and qj{n), which we 
denoted by Dp. Partitions of two vertical rectangular cells were considered for 
various cell widths. It was found that agreement was best when all cells were 
similar, but overall agreement was generally good {Dp < 0.15). In particular, 
uniform rectangular partitions gave exact agreement {Dp = 0). This was un- 
derstood to be a consequence of the fact that such partitions are in fact Markov 
partitions for the baker map. In general, however, the agreement between the 
exact and approximating processes will depend upon how strongly the system 
is mixing. Since the baker map is not only mixing but also a Bernoulli system, 
it is not clear whether such good agreement is typical of most mixing systems. 

While the time evolution of a coarse grained observable under a deterministic 
map may not be Markovian, for mixing maps this evolution is at least approx- 
imately Markovian in the sense that both {Yn} and {Yn} have the same dis- 
tribution initially and in the distant future. (If the map is invertible then they 
also have the same distribution in the distant past.) Thus, the two processes 
do not differ qualitatively and in some cases agree quite well quantitatively. 

We have not addressed the problem of obtaining closed form deterministic laws 
for macroscopic variables. This may be possible for "many to one" observ- 
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ables corresponding to strongly nonuniform partitions. In general, however, 
such macroscopic laws are expected to be only approximately deterministic, 
as befits any realistic model. 
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Appendix 



Here we prove the claim that every uniform rectangular partition is a Markov 
partition for the baker map. Note that this is not a true partition since the 
boundaries of the cells overlap; however, the boundaries are of measure zero 
and may easily be modified to form a true partition. 



PROOF. (Theorem 4) The collection {Rij} of rectangles will be a Markov 
partition provided (x, y) G int Rij and 0(x, y) G int Ri'j' imply that (f>(W^. , {x, y)) 3 
WR^,^^X<Pi.x,y)) and 0(l^R.^.(a;,|/)) C l^^^,^,(0(a;,|/)), where W^J|^ ^, (x, y) and 
Wji. .{x,y) are the intersections with Rij of the unstable and stable mani- 
folds, respectively, at {x,y) [14]. For the baker map we have W^, X^iU) = 

[^^Ji] X {y} and Wl^^^^{x,y) = {x} x [^,^]. We consider the following 
three cases separately: 

Case 1: i^ < I 

Since a; < -^ < |, {x',y') = (j){x,y) = (2x, |) and so 



<^K,,(a:,l/)) = [^,i]x{f} 
<l>{W^Jx,y)) = {2x}x[^,^] 



and 



WLM{x,y)) = {2x}x[^^,^]. 



i-l ^ ^. ^ i 1 



We now must determine i',j' in terms of i,j. Note that if Vj- < x < -^ — 



M ^ •" ^ M 2M 



then ^^ < x' < \f-, so i' = 2i - 1. Similarly, ifT7-5T7<a;<-^ then 
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i' = 2i. Note also that ^ < y < j^ implies ^^ ^ v' ^ 4n' ^^ ^' ~ I ^^ i i^ 



even and j' = ^ when j is odd. Thus, 



^<ig-W ■ Wl(<l>ix,y)) = [^,'-^]x{y^} 



2J ' 



J even : W^ M{^,y)) = {2a;} x [i^, ^] , 
J odd : W-l (0(a;,y)) = {24x[i^,i^]. 



Case 2: ^ < 1 < -^ 



For this case M is necessarily odd and i = ^^^. We consider the cases x < ^ 
and X > -^ separately. For both cases we have 

<p{W'nJx^y)) = [l~i^,l]x{^} U [0,^]x{^}. 



Case 2a: x < | 

Since x < ^, {x',y') = (f>{x,y) = {2x, |). From the fact that ^^ < x < | it 
follows that ^ < x' < 1, so i' = M and 

Wl,M{x,y)) = [M^,l]x{y^}. 



For the stable manifold we have [W^. .{x,y)] = {2x} x [^^j ^], so 



J even : W^4,^,,(0(x,i/)) = {2x} x [y, ^fe] , 
J odd : W^^,J<P{x,y)) = {2x}x[^,^], 

since j' = | when j is even and j' = ^^ otherwise. 

Case 2b: i < X 

Since x > |, (x', y') = 0(x, y) = (2x— 1, ^^). From the fact that | < x < 
it follows that < x' < -^, so i' = 1 and 



For the stable manifold we have [W^^ . (x, y)) = {2x — 1} x [^ + I' 27v + il • 
From the fact that ■^ < y < ^^ it follows that ^ + |<y'<2iv"'"i- Thus, 
f = ^^^ ii j + N is even, and j' = ''^ o^^ ii j + N is odd. This gives us 



M+l 
2M 
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j + iVeven : 1^^ (0(x,y)) = {2a: - 1} x [^ + i, ^ + i] , 
J + AT odd : W%,U{x,y)) = {2x-l}x[i^ + \,i^ + \]. 



Case 3: I < ^ 



Since 2: > jg > |, (a;', y') = (j){x, y) = {2x — 1, ^ji) and so 



i{wi,Jx,y)) = {2x-l}x[^ + \,^ + \] 
and 

1^4 .,(0(x,i/)) = {2a:-l}x[^,|:]. 

Note that if ^ < a; < -^-^^then^-l < x' < ^-1, so i' = 22-1-M. 
Similarly, '^^ ^ - ^ < x < ^ then i' = 2i - M. Note also that j' = ^ if 
j + iV is even and j' = •^'"'"•^"'"^ when j + iV is odd. Thus, 

j + iVeven : W^^ ,(0(x,y)) = {2a; - 1} x [^ + i, ^ + i] , 
j + iVodd : 1^4 ,(0(a;,y)) = {2a;-l}x[i^ + i,^ + i]. 



This completes the proof. □ 
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